МІНІСТЕРСТВО ОСВІТИ І НАУКИ,МОЛОДІ ТА СПОРТУ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ “ЛЬВІВСЬКА ПОЛІТЕХНІКА”
ІКТА
кафедра ЗІ
З В І Т
до лабораторної роботи №1
з курсу: «Цифрова обробка сигналів»
на тему: «Дискретне перетворення фур’є та його застосування для спектрального аналізу сигналів»
Варіант №6
Мета роботи – ознайомлення із математичним апаратом опису сигналів у частотній області, змістом дискретного перетворення Фур’є та його застосуванням для спектрального аналізу реальних сигналів.
ЗАВДАННЯ
1. Навести аналітичний вираз та обчислити спектральні коефіцієнти періодичного сигналу, одержаного шляхом двонапівперіодного випрямлення гармонічного коливання, із параметрами табл.1 . Добрати параметри ДПФ для спектрального аналізу періодичного сигналу, щоб забезпечити вимоги в табл. 1. Показати графіки часової функції сигналу і його спектра.
Таблиця 1
№ п/п
АмплітудаAm, В.
Період коливання T0, с.
К-стьспектральних
коефіцієнтів
Роздільчаздатністьпо частоті ΔF, Гц
6
2
0,6
6
1/6
clc; % очистити командне вікно
clear all; % звільнити пам'ять робочого середовища
%%%%%%%%%% ВАРІАТИВНА ЧАСТИНА %%%%%%%%%%%%%%%%%%%%%%
Am=2;% Амплітуда гармонійного коливання
k=6; %кількість спектральних коефіцієнтів
Tk=0.6; % період гармонійного коливання
dF=1/6; % роздільча здатність по частоті
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
T0=Tk/2; % період сигналу одержаного шляхом двонапівперіодного випрямлення
Fmax=k*1/T0;% k-та гармоніка у спектрі сигналу (максимальна частота)
Ts=1/(2*Fmax); % період дискретизації за теоремою Котельникова
T=1/dF; % інтервал спостереження ДПФ, що забезпечує роздільчу здатність dF
t=0:Ts:T-Ts; % N дискретних моментів часу на інтервалі спостереження
s=abs(Am*sin(2*pi*t/Tk));% N значень сигналу у дискретні моменти часу
y=fft(s); % швидке ДПФ,
% y приймає N дискретних спектральних значень у діапазоні [0..Fs]
yy=fftshift(y); % перестановка правої і лівої частин для відображення у діапазоні [-Fs/2..Fs/2]
figure(1);
tau=0:Ts/8:T-Ts;
plot (tau,abs(Am*sin(2*pi*tau/Tk))); axis([0 Tk 0 Am]); hold on
stem(t,s); hold off % часове представлення сигналу
figure (2);
xx=abs(yy/(length(yy))); % амплітудний спектр сигналу (з множником 1/N)
f=-1/(2*Ts):dF:1/(2*Ts)-dF; % частоти, що відповідають N дискретним спектральним значенням
stem (f,xx); % частотне представлення сигналу
figure(3);
s=[s s s];
plot (-t(length(t))-Ts:Ts:2*t(length(t))+Ts, s);axis([-Tk 2*Tk 0 Am]); hold on
stem (-t(length(t))-Ts:Ts:2*t(length(t))+Ts, s); hold off;
Результати роботи програми:
/ /
Рис.1 Часова функція періодичного сигналу Рис.2 Спектр періодичного сигналу
/
2. Навести аналітичний вираз спектральної густини експоненціального імпульсуs (t)=Am×exp(-|a×t|), параметри якого наведено в табл. 2. Добрати параметри ДПФ для спектрального аналізу імпульсного сигналу, щоб забезпечити вимоги даного варіанту. Показати графіки часової функції сигналу і його спектра.
Таблиця 2
№ п/п
Амплітуда Am, В.
Сталазгасання a, с-1
Частотнийінтервал, Гц
Роздільча здатність по
частоті ΔF, Гц
6
2
0,6
1,5
0,05
clc; % очистити командне вікно
clear all; % звільнити пам'ять робочого середовища
%%%%%%%%%% ВАРІАТИВНА ЧАСТИНА %%%%%%%%%%%%%%%%%%%%%%
Am=2; % Амплітуда імпульсу
a=0.6; % стала згасання експоненційного імпульсу
Fmax=1.5; % частотний інтервал
dF=0.05; % роздільча здатність по частоті
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ts=1/(2*Fmax); % період дискретизації за теоремою Котельникова
T=1/dF; % інтервал спостереження ДПФ, що забезпечує роздільчу здатність dF
t=0:Ts:(T-Ts); % N дискретних моментів часу на інтервалі спостереження
s=Am*exp(-abs(a*t)); % N значень сигналу у дискретні моменти часу
y=fft (s); % швидке ДПФ,
% y приймає N дискретних спектральних значень у діапазоні [0..Fs]
yy=fftshift (y); % перестановка правої і лівої частин для відображення у діапазоні [-Fs/2..Fs/2]
figure(1);
tau=0:Ts/8:T-Ts;
plot (tau,Am*exp(-abs(a*tau))); hold on
stem(t,s); hold off ...